Multivalent ion effects on electrostatic stability of virus-like nano-shells 
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Electrostatic properties and stability of virus- like nano-shells are examined in mixed ionic solutions 
with monovalent and multivalent ions. A theoretical model based on a thin charged sphere and 
multivalent ions within the "dressed multivalent ion" approximation is compared with extensive 
implicit Monte-Carlo simulations and shown to be accurate for positive or low negative surface charge 
density of the shell, and for sufficiently high (low) monovalent (multivalent) salt concentrations. 
Deviations between simulation and theory occur due to many-body effects leading to formation of 
a strongly correlated layer of vicinal ions. Phase diagrams of electrostatic (osmotic) pressure show 
positive and negative pressure values, corresponding to outward and inward facing forces on the 
shell, providing an explanation for high sensitivity of viral shell stability on the ionic environment. 
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Electrostatic interactions contribute a significant part 
to the total free energy of charged viral capsids and/or 
virus-like shells In this context multivalent ions be- 
have rather differently from monovalent ones. Some of 
the major functions of multivalent metal ions such as 
Mg 2+ , Zn 2+ , and Ca 2+ are surmised to enhance the ca 
sid stability, enable their conformational alterations 
or assist in the genome packaging 0], while multiva- 
lent polyamines like spermine (Spm 4 ~M are required for 
a proper stabilization of viruses d Hf. Though Ross- 
man and Eriksson, already a while ago, clearly associated 
swelling of capsids with the removal of multivalent ions 
the exact mechanism of electrostatic capsid destabi- 
lization remains to a large extent elusive. Here we sug- 
gest a rather universal mechanism that could underline 
the action of multivalent ions on viral capsid stability. 

In order to properly understand the nature of elec- 
trostatic interactions in viruses, we investigate the ef- 
fects of highly charged bathing solution ions on the cap- 
sid stability via their spatial distribution and the as- 
sociated (osmotic) pressure. To accomplish this goal 
and analyze the distribution of highly charged ions 
in the vicinity of charged capsid shells, the Poisson- 
Boltzmann (PB) paradigm, valid only for weakly charged 
or highly screened Coulomb systems, has to be prop- 
erly amended 0- This is achieved in a twofold manner: 
by formulating an analytical "dressed multivalent-ion ap- 
proach" as well as by combining it with exhaustive 
Monte-Carlo (MC) simulations within a simple model in 
order to properly capture the salient characteristics of 
electrostatic effects on capsid stability. 

Consider a closed spherical shell of radius R and sur- 
face charge density a, Fig. QJi, which may serve as a rea- 
sonable model for a wide range of viral capsids such as 
brome mosaic virus, cowpea chlorotic mottle virus, etc., 
[H, as well as synthetic nano-shells 11, 12|. whose 



effective radius and surface charge density can both be 
extracted from experimental data. For typical viruses, 
these values fall in the range 2.5 nm < R < 30 nm with 
thicknesses between 2 and 4 nm and |er| < 0.4 eo/nm 2 . 
There is no clear-cut preference for inner-outer capsid 
surface charge densit y si gn segregation and in fact both 
signs are admissible [lOj. The shell is considered to be 
permeable to ions and immersed in an ionic solution com- 
prising a monovalent 1 : 1 salt of bulk concentration no 
and a multivalent q : 1 salt of bulk concentration Co, 
where q is the charge valency of multivalent ions. 

The above model presents an asymmetric multicom- 
ponent Coulomb system which is not amenable to exact 
analytical treatment. We therefore resort to the "dressed 
multivalent-ion approximation" Q, which can be de- 
rived as an approximate theory from a systematic field- 
theoretical formalism in situations where the ionic mix- 
ture is highly asymmetric, i.e., q > 1. In this approach, 
the degrees of freedom associated with monovalent ions 
are integrated out from the partition function, result- 
ing in an effective formalism incorporating only screened 
Debye-Hiickel (DH) interactions between the remaining 
"dressed" multivalent ions and fixed (macroion) surface 
charges On this level, the problem is still a many- 
body one but it can be treated with fast and efficient MC 
simulations as one needs to include only the dressed ions' 
degrees of freedom. In order to proceed with an analyti- 
cal theory, however, it turns out that the grand-canonical 
partition function of the system can be expanded in 
terms of the fugacity (concentration) of multivalent ions, 
which is indeed small in most cases of experimental in- 
terest 0, To the leading order, this procedure then 
leads to the dressed multivalent-ion theory containing the 
lowest order contribution from the interaction between 
dressed multivalent ions and fixed macroion charges in 
the system. The dressed multivalent-ion theory is thus 
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FIG. 1: (Color online) a) Schematic view of a spherical shell of radius R and surface charge density a immersed in a mixed 
solution comprising mono- and multivalent salt ions (multivalent ions are shown by bigger red spheres), b) Rescaled density of 
multivalent ions, c(r)/co, in the vicinity of a charged shell with radius R = 10 nm as a function of the radial distance r from the 
center of the shell. Symbols show MC simulation data and lines the analytical dressed multivalent-ion predictions for a = +0.4 
eo/nm 2 , q = 4, and various k as indicated on the graph, c) Same as (b) but for a — —0.05 eo/nm 2 , and k and q values as 
shown on the graph, and d) for k = 1 nm -1 , q = 4, and different values of (negative) a. In all cases, we have qno = 4 mM. 



related to the strong- coupling theory, introduced origi- 
nally for salt-free systems containing only multivalent 
counterions Q. It has been tested extensively against 
all-ion (explicit) MC simulations and has been shown to 
be a valid theory for moderately large values of q and 
monovalent salt concentration no , and for sufficiently low 
concentration of multivalent salt cq [8, 14 1. 

For a charged shell, one can express the dressed 
multivalent-ion free energy of the system as follows [l,[l4[ 

l3T = l3T DH -c a [ dre-^o^W +0(c 2 ); (1) 



where j3 = l/(fcsT) and the first term is the self-energy 
of the charged shell in the DH limit, i.e., in the absence 
of multivalent ions, which can be calculated as in Ref. [15. 
The integral in the second term corresponds to the par- 
tition function for the single-particle DH interaction of 
a multivalent ion with the charged shell, with <y5_D.fr (r) 
being the screened electrostatic potential of the shell, 
which follows simply from the DH equation V 2 (fiDH (r) = 
k 2 (Pdh(j) solved in the presence of the charged perme- 
able shell. Note that the inverse Debye screening length 
k depends on the bulk concentrations of both mono- 
and multivalent ions as k 2 = 87t^b("-o + <Z c o/2), where 

= /3eo/(47T££o) is the Bjerrum length. 

The (number) density of multivalent ions is then 



c(r) = c exp ( - f3qe ip DH (r)) , 



(2) 



and the electrostatic (osmotic) pressure acting on the 
shell is obtained as P = —dJ-/dV, or 



dV 



(3) 

where V = 47ri? 3 /3 is the shell volume, and the partial 
derivatives are taken at fixed total charge Q a = AiraR 2 . 
The pressure can thus be divided into two parts, one 



stemming from the sclf-cncrgy of the charged shell Pdh, 
and the other one from the multivalent ions. The for- 
mer contribution is purely positive but the latter can be- 
come negative as we explain later. The explicit forms of 
<Pdh(t) = (fDnir) and Fdh can be found in Ref. [15. 

In order to assess the validity of the dressed 
multivalent-ion theory detailed above we have used ex- 
tensive MC simulations in which we simulate a large 
collection of mobile (dressed) multivalent ions near a 
charged permeable shell with the monovalent ions be- 
ing treated implicitly. This means that multivalent ions 
as well as the fixed charges on the shell interact through 
the screened DH interaction pair potential. The monova- 
lent salt concentration, no, is varied in the range between 
30 mM to 300 mM and we typically take a small multi- 
valent salt concentration Co of a few mM in accordance 
with virus experiments 0, EH- We use canonical MC 
simulations for dressed multivalent ions (assumed to be 
confined in a sufficiently large outer cubic volume with 
the charged shell located at the center), but employ an it- 
erative method to produce a constant equilibrium "bulk" 
concentration of multivalent ions of a given value at large 
separations from the capsid [jq . which thus effectively 
simulates an open (grand-canonical) system. In what 
follows, we shall take positively charged multivalent ions 
with q = 2, 3, 4, the ambient temperature T = 300 K, 
and the dielectric constant as that of water, e = 80. 

In Fig. [TJ we show the rescaled number density profile 
of multivalent ions in the vicinity of a charged shell of ra- 
dius R = 10 nm as a function of the radial distance from 
the center of the shell. In certain regimes of parameters 
the analytical prediction from the dressed multivalent- 
ion theory [Eq. ([2]), shown by solid lines] performs very 
well and closely matches the MC simulation results (sym- 
bols). For positively charged shells, a > 0, there is always 
an excellent agreement showing a drop in the density of 
ions close to the shell surface (Fig. [TJd). For negative <r, 
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however, multivalent ions are attracted and accumulate 
more strongly at the shell surface. There is still a wide 
range of parameters where theory and MC simulations 
quantitatively agree (Figs. QJ and d, e.g., for divalent 
ions and the DH screening length k = 0.58 nm -1 corre- 
sponding to a 30 mM monovalent salt concentration) . In 
fact, as shown in Fig. by increasing q, the multiva- 
lent ion density at the surface increases but at the same 
time deviations between the simulations and the analyt- 
ical predictions tend to grow. This implies an incipient 
role for multi-particle correlations which are absent in the 
theory. In all cases, where there is a deviation, the theory 
appears to overestimate the density of ions in the vicinity 
of the shell. For large q, an improved agreement can be 
achieved by increasing the screening parameter which re- 
duces multi-particle interactions; e.g., for tetravalent ions 
the analytical prediction works very well for k = 1 nm -1 , 
corresponding to 100 mM of 1:1 salt (black circles). 

A similar trend is observed when the sign of surface 
charge density of the shell, a, turns negative. In this 
large increase is obtained (by up to a few orders 
of magnitude relative to the bulk concentration, Co) in 
the surface density of tetravalent ions when the surface 
charge density is only modestly increased from a = —0.05 
eo/nm 2 to a = —0.4 eo/nm 2 (Fig. [Hi). The theory again 
works at sufficiently low \a\ and breaks down at larger 
values. Note that the deviation between the theory and 
simulations can be tuned further by changing the screen- 
ing parameter: the validity of the theory for a = —0.4 
eo/nm 2 is restored (not shown) [l6| if the screening pa- 
rameter is increased to K = 1.73 nm , equivalent to 
300 mM of 1:1 salt. The break-down of the theory is 
more dramatically observed at n — 0.58 nm -1 with both 
limits falling well within the experimentally relevant val- 
ues. Thus, the applicability of the above analytical pre- 
dictions depends on all three parameters, q, a and k [l6(. 

While the density of the multivalent ions across the 
shell surface can be understood within the phenomenol- 
ogy of the dressed multivalent-ion theory, the behavior 
in the lateral directions along the shell surface can be 
less so. The pair distribution function (pdf) of multi- 
valent ions along the surface (Fig. [2^) indeed shows a 
noticeable correlation hole when plotted in terms of the 
lateral distance between multivalent ions, which results 
from the repulsion between ions residing at the vicinity 
of the shell. The size of the correlation hole, a±, is con- 
sistent with the naive estimate a± ~ ^Jqea/\a\ (being 
around 2.2 nm for the parameters in Fig. (2^), which can 
follow from a local electroneutrality argument for surface 
ions Q- More remarkably, however, in the regime where 
the accumulation of the ions at the shell surface is highly 
enhanced and the theory breaks down (Fig. [T]), we find a 
short-range ordering of multivalent ions along the surface 
clearly exhibited by a pronounced correlation peak in the 
pdf of multivalent ions. This behavior illustrates the fact 
that the failure of the dressed multivalent-ion theory is 
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FIG. 2: (Color online) a) Normalized pair distribution func- 
tion of multivalent ions from MC simulations for R = 2.5 nm, 
a = —0.8 eo/nm 2 , k = lnm -1 , q = 4. b) A typical simulated 
snapshot showing multivalent ions in the vicinity (within lA) 
of the shell for R = 2.5 nm, a — —1.6 eo/nm 2 , k = 1.73 nm -1 , 
q = 4. For these parameters, the shell experiences a total neg- 
ative pressure of (a) —0.28 feT/nm 3 and (b) —1.31 teT/nm 3 . 
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FIG. 3: (Color online) a) Rescaled total electrostatic pressure 
Pr = sgn(P') x log 10 (1 + 100 |P'|) (where P' = P/P with 
Po — 1 feT/nm 3 ) from simulations (symbols) and the dressed 
multivalent-ion theory (lines) as a function of a for R — 2.5 
nm, q = 4, Co = 1 mM, and different values of n. b) Same 
as (a) but here we show different components of the total 
pressure due to shell's self-energy, Pdh (green, upper lines) 
and multivalent ions, P q (red symbols and red, lower lines). 

driven by the formation of a strongly correlated surface 
layer of multivalent ions, which incorporates many-body 
effects and is thus absent from the presented theory. 

The electrostatic pressure acting on the shell also 
shows interesting variation depending on the system pa- 
rameters. In Fig. [2k we plot the rescaled pressure, Pr, 
as a function of the surface charge density of the shell 
for tetravalent ions and R — 2.5 nm (the rescaling is 
done in such a way that the negative and positive val- 
ues of pressure can be shown in a log-linear scale and 
points of zero pressure are correctly represented). For 
positively charged capsids, the pressure is always posi- 
tive. In this case there is excellent agreement between 
analytical prediction (blue lines) and simulations. For 
negatively charged capsids, both the theory and simu- 
lations turn out to give negative pressure for screening 
parameters in the vicinity of the physiological value or 
smaller, indicating that the shell is being compressed. In 
this case the theory again tends to overestimate the mag- 
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FIG. 4: (Color online) "Phase diagrams" showing regions of 
positive and negative electrostatic pressure acting on the shell 
in terms of the system parameters, q and a (shown in linear 
scale) as well as R and n (shown in logarithmic scale). Solid 
lines show the line of zero pressure. In a) we have fixed R — 
2.5 nm and k = 0.5811m" 1 , in b) we have fixed q — 4 and 
k — 0.58 nm - , in c) q — 4 and R = 2.5 nm, and in d) we have 
fixed q = 4 and the total shell charge Q a — AivaR 2 = —126 eo 
(e.g., we have a ~ —0.1 eo/nm 2 for R = 10 nm). 

nitude of the negative pressure and shows deviations from 
the simulations (compare symbols and lines in Fig. [3jt) . 
When the contribution due to the self-energy of the shell 
Pdh is subtracted, one obtains the pressure imparted 
on the shell by multivalent ions, P q , which is shown by 
red symbols (simulations) and lines (theory) in Fig. [5Jd. 
Clearly, the negative pressure thus results purely from 
multivalent ions and then only when the shell is overall 
negatively charged. This negative pressure is directly re- 
lated to the attractive interactions mediated by highly 
charged ions between equally charged surfaces [1J|; it in 
fact results from summing up attractive forces along the 
outer and inner surfaces of the charged shell that have 
a component of the force directed towards the inside of 
the shell [r|. The sign of osmotic pressure acting on 
the shell is thus straightforwardly related to correlation 
attractions as seen in other strongly charged systems Q ■ 

Note that an important factor determining the sign of 
the pressure is the bulk concentrations of multivalent as 
well as monovalent ions, as codified by k, which can lead 
to a change of sign of the total pressure due to an increase 
in the negative pressure from the multivalent ions. In 
Fig. 0] we thus plot "phase diagrams" of the pressure 
acting on the shell as obtained from our MC simulations 
in terms of the parameters q, a, n, and R. We show 
regions where the pressure is negative and/or positive, 
with the black solid line indicating the separatrix of zero 
pressure. It is thus clear that by changing the system 
parameters one can reach a point where the electrostatic 
self-rcpulsion of the capsid is overcome and gives way to 



an inward force isotropically compressing the capsid. 

The sign and magnitude of the pressure are both cru- 
cial in assessing the capsid stability towards shrinking 
and/or swelling. Sufficiently large negative pressure can 
lead to inward buckling of the capsid through a sequence 
of shrinking clastic instabilities |17| . expected to occur 



when negative pressure is large enough, exceeding a crit- 
ical value of P c ~ —5 [atm] x (30 [nm]/i?) 3 in magni- 
tude. Our results indicate that in the range of a ~ —0.5 
eo/nm 2 and Co ~ 1 mM, this type of destabilization and 
inward buckling may take place for large shells with ra- 
dius R > 150 nm, an estimate that could however be 
lowered significantly by increasing the concentration of 
multivalent ions [ly]. On the other hand, the change of 
sign of the pressure on complete removal of multivalent 
ions that turns the negative pressure to positive, would 
lead to swelling of the capsid and possibly outward buck- 
ling, destabilizing it towards swelling |6(. Multivalent 
ions can thus be implicated in both types of instabilities, 
leading to larger effects for larger sized capsids. Inter- 
estingly in this context, larger viruses such as mimivirus, 
megavirus, etc., with sufficiently large radii all seem to 
have either an external or an internal lipid membrane 
[l8| that could serve as a local reservoir for multivalent 
ions acting in cither direction. In general, local multiva- 
lent salt concentration gradients generated by a semiper- 
meable membrane, could also play an important role in 
regulating the stability of other enveloped viruses. 

Negative osmotic pressure also leads to important con- 
sequences for self-assembly of viral shells. Clearly for 
P > 0, the size of the capsid spontaneously increases, 
precluding a stable state with a finite R. On the other 
hand for P < 0, there could be two conceivable out- 
comes. In this regime strong electrostatic coupling will 
provide effective attractive interactions between equally 
charged capsomcrs, thus promoting their aggregation but 
not necessarily also their self-assembly into a fully formed 
capsid. A disordered aggregate might be a more stable 
configuration then an ordered spherical shell, indicating 
that salt and multivalent ion concentration could affect 
the strength of effective lateral interaction between cap- 
somers and thus interfere with self-assembly. There are 
indications, mostly for divalent salts, that this is indeed 
what happens and that solution conditions (no of salt, q, 
and Co of multivalent ions) do play an important role in 
the assembly /disassembly process 
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